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Abstract. A model is suggested for the bursting x-ray 
pulsar, GRO J1744-28, where the compact object of the 
pulsar is proposed to be a strange matter star with crusts. 
It is suggested that two envelope crusts shield the two 
polar caps of an accreting strange stars owing to strong 
polar-cap magnetic fields. GRO J1744-28 might not have 
a unified crust, but two polar crusts. Bursts are the result 
of phase transition of polar crusts if one accreting crust 
is heavier than that the Coulomb force can support. It 
is shown that there should be phase-lags when bursting 
crust is expanded about one hundred meters. Many calcu- 
lated properties of this model, such as the phase-lags, the 
bursting luminosity and the characteristic photon energy, 
are in agreement with observations. 
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1. Introduction 

GRO J1744-28 is the first known X-ray source both burst- 
ing and pulsating (see Rappaport & Joss 1997, Aptekar et 
al. 1998, and references therein). The source is in a bi- 
nary system, the X-ray pulsation period is 467 ms, the 
orbital period is 11.8 days, the mass function of this sys- 
tem is 1.31 x 10~ 4 M Q , the time interval between bursts is 
about 30 minutes, and the duration of the bursts is about 
5 seconds. There is not any distinct difference between the 
spectra of the bursts and persistent emission, and the typ- 
ical photon energy is about 14 keV. When bursting, there 
are post-burst dips in the persistent flux and phase-lags. 



The phase of the pulse profile during bursts lags the phase 
prior to the bursts by 70°-50°. 

Any theoretical model should make clear that, such 
as, how matter accretes onto the surface of GRO J 1744- 
28, how the bursts produce, why there are phase-lags 
and post-burst dips, etc. Usually, people believe that the 
source could be well explained by a complex accreting pro- 
cess onto neutron star. One of this kind of model believe 
that the bursts are the result of accretion instability (Type 
II X-ray bursts, Lewin et al. 1996), which has been stud- 
ied in the Rapid burster (MXB 1730-335). On the oher 
hand, Bildsten & Brown (1997) show that the bursts ob- 
served during the peak of the outburst are most likely not 
the thermonuclear origin, i.e., GRO J1744-28 could not be 
Type I X-ray burster. 

However, Cheng et al. (1998) have shown the possi- 
bility that the burst energy is the result of phase transi- 
tion of accreted matter from normal hadron to three-flavor 
strange quark matter. Hence, they concluded that GRO 
J1744-28 might be an accreting strange star. In their dis- 
cussion, they address that GRO J1744-28 is a strange star 
with a normal matter crust. When the accreted mass ex- 
ceeds some critical mass, the crust breaks, resulting in the 
phase transition and the energy release. 

Xu & Qiao (1998) argued that the normal matter crust 
is not necessary for the strange star model of pulsars. If 
a strange star (Alcock et al. 1986) forms soon after a su- 
pernova, a magnetosphere mainly composed by would 
soon be established, and the strange star could act as a 
radio pulsar. There are some advantages if radio pulsars 
are bare strange stars rather than neutron stars or strange 
stars with crust. While strange stars are in binaries, there 
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could be two envelope crusts shielding the two polar caps 
of the accreting strange stars. 

Following this idea, we propose a model for GRO 
J1744-28, which is similar to that by Cheng et al.(1998). 
We argue that GRO J 1744-28 might not own a unified 
crust, but two polar crusts. Bursts are the result of phase 
transition of polar crusts if one accreting crust is heav- 
ier than that the Coulomb force can support. This phase 
transition process could be unstable, similar to the case 
discussed by Horvath & Benvenuto (1988). It is shown 
that there should be phase-lags when bursting crust is ex- 
panded about one hundred meters, a scale length in order 
of the crust thickness. 

About the phase-lag, Miller (1996) has suggested that 
when accretion rate rises, the flattened accretion flow shift 
to a different set of field lines. Therefore, the time of pulse 
maximum would change accordingly, and a phase-lag re- 
sults. Our model is different from this scenario. 

2. General Description of the Crust over Strange 
Quark Matter 

As s quark is slightly heavier than u and d quarks, there 
are u, d, s quarks and electrons for the chemical and elec- 
trical equilibrium of strange quark matter. Since the quark 
matter are bound up through strong interaction, and the 
electrons are held to the quark matter electrically, the elec- 
trons' distribution would extend beyond the quark matter 
surface several hundred fm, and there is a strong elec- 
tric field near the quark surface. Using the Thomas-Fermi 
model, one can obtain the electric field E as a function of 
the height above the quark surface z (Xu & Qiao, 1998, 
cqu. A5-A6), 



E = - 



dV 

dz 



9V„ 



( 

1.42X10 21 



Vg-z+4) 2 

Vm-\ 



(1) 



where, the direction of the electric field is outward, and 
the potential energy of quark matter V q has been chosen 
to be 20 MeV, z n = z/(Kr n cm). 

As matter accreting on to the polar caps, the strange 
quark matter should repulse atomic nuclei, and two crusts 
are formed in the two polar caps. Assuming the density of 
the crust is /?cmst, the distance I between nuclei is 



I = ( 



Axm p -.1/3 

Pcrust 



where m p is proton mass, A is the mass number of nucleus, 
A56 = A/56, pu = Pcmst/lO 1 ^ cm" 3 . If p crust approxi- 
mates neutron drip density, p\\ ~ 1 (Glcndcnning et al. 
1995), then I is order of the electric gap height (Alcock et 
al. 1986). 

By simulation, Huang & Lu (1997) argued that the 
electric gap between crust and strange quark matter is too 
small to support the crust before the density at the bottom 
rises to the neutron drip density. Because of the repulsion 
between nuclei and strange quark matter, there should be 
a charge separation near the bottom of crust. Hence there 
would be an effective positively charged column density 
cr c , which could be estimated as 

= V ■ fr 

~ 2.7 x 10 23 V p 2 1 { 3 A^ /3 Z 26 e, 

where 77 ~ 1 is a coefficient denoting the effective (posi- 
tive) chargement, Z is charge number of atomic nucleus, 
Z 2 6 = Z/26, e is the charge of electron. Hence, by equilib- 
rium between electrical and gravitational forces, the col- 
umn density cr m of crust can be written as 

._ a c ER 2 
° m ~ GM, 

(1 2211+4)^ 11 A 56 Z 2S M 1 H 6 g'Cm 

- ^^VPTA-^Z 2 ,M-^RI M .cm-2, 

(2) 

where M* and R are the mass and radius of strange quark 
matter, respectively, Mi — M»/(1M Q ), G is the gravita- 
tional constants. The height H crust of the crust is 



1.1 x 10 4 r ? p 11 1/3 A 56 2/3 Z 26 M 1 - 1 Rl cm (for z u = 



For a spherical crust, the total mass M crU st should be 
M crust = 47ri? 2 cr m 

29xl0~ 5 2/3 4-2/3 7 , r l R 4 M 
(1.2^11+4)^ 11 Ai6 Z^M X R 6 Mq, 

which is a typical value of simulated results (e.g. Huang 
& Lu 1997). For zu — (the length scale of strong inter- 
action is ~ zn = 10~ 2 ), we get the maximum values of 
column density er max and of spherical total mass M max , 



<W = 1.4 x 10- 1 »r 1 p 2 1 / 1 3 A- 2/3 Z 26 M- 1 R' 6 M • cm 



lr>2 



56 

M max - 1.8 x lO- 5 r]p 2 1 { 3 A-* /3 Z 26 Mi 1 Rt Mq. 

(3) 

If the mass of the crust is known, the distance be- 
tween the crust and strange quark matter should be (from 
equ.(2)) 

An 1/2 -1/2 1/3 4-1/3 ,7-1/2 7.^-1/2 D OD 

z lx - 4.(V /Z cr 18 pl x A 56 < Z^ M 1 Re -3.3, 



1.0 x 10~ u A$ Pl{ 



1/3 



cm, 
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where ais = er m /(lO~ 18 M0cm~ 2 ). It is shown from above 
equation that the crust becomes heavier as matter ac- 
cretes, the thickness of the electric gap becomes smaller. 

3. Bursting Process above the Polar Cap 

As matter accreting to the polar cap crusts, the column 
density of crusts could be more and more higher, until a 
catastrophe of crust takes place when cr m > er max , because 
the strong interaction length is of order 1 fm. This crush 
of crust could trigger a burst process in GRO J1744-18. 

When a crush happens, the strange quark core should 
swallow part of normal matter at the bottom of the crust. 
This process could not be unstable (possibly in a detona- 
tion mode), which might be similar to that discussed by 
Horvath & Benvenuto (1988). 

Before going to some details of the model for GRO 
J1744-28, we would like to briefly review the three types 
of particular bursts. There are two types of X-ray bursts, 
which we call them as type I and type II bursts (Lewin et 
al. 1993). The ratios of steady emission to time averaged 
burst emission in type I bursts, a\, and in type II bursts, 
an, are 

10 20 crgs-s~ 1 



ai 
an 



10 ls or§s-s _1 

1.3xl(r"crgs-s^ 

1.3xlO ; "crgs-s- 1 



100, 

1, 



respectively, while, the Eddington luminosity L-&dd ~ 
10 38 ergs • s _1 . For the bursts in GRO J1744-18, in the 
early 1995 Dec, when the burst intervals were ~ 3 min- 
utes, the 3ct upper limit of the ratio of integrated energy 
in the persistent flux to that in bursts was ~ 4 (Lewin et 
al. 1996), i.e., 

"1744 < 4. 

There are three main energy sources in the X-ray emis- 
sion by unit nucleon. First, the gravitational accretion en- 
ergy, eq ■ Second, the thermonuclear reaction energy from 
pure hydrogen into helium, en, and from helium into iron- 
peak elements, £hc- Third, the energy liberated in the 
phase transition from normal hadron into strange quark 
matter, e s . They are 

e G ~ 140MeV/nucleon, 

en ~ 6.7MeV/nucleon, 

£Hc ~ 1.7McV/nucleon, 

e s ~ 30MeV/nucleon. 

For type I X-ray bursts caused by nuclear flash, the 

^ — to ~ For type 



observed ratios ct\ arc from 



II X-ray bursts, which are the results of accretion instabil- 
ity, the observed ratios an ~ ~ 1 if the time averaged 
accretion rates during bursts and between bursts are sim- 
ilar. 

While, for the busting X-ray pulsar, GRO J1744-28, 
there should be a phase transition from hadron to three- 
flavor quark matter when the crust crashes. Nevertheless, 
the energy release from the phase transition could trigger 
the thermonuclear reactions in the crust. Hence, the ratio 
of steady emission to time averaged burst emission could 



be 



3.6, which is comparable to the observed 



«1744- 

The polar cap area of the crust, A p , is 



A p = 



nR 3 



where, r c = min[r c , r m ], and r c is the radius of light cylin- 
der, r m is the radius of magnetosphere (Rappaport & Joss, 
1997), 



cP 
2tt 



2.2 x 10 9 cm, 

4/7 p 12/7 



= 2.8 x 10 7 BifR, 



M\ 



-1/7 



-2/7 



cm. 



Where B is the magnetic field, B\\ = _B/(10 n gauss), M 
is the accretion rate, Ms— ir ,_ 8 Jf — =t, P = 0.467s is the 
period of rotation, c is the velocity of light. Hence (let 



A p = 1.1 x lO L "B~^'R t 



< I „-4/7 p 9/7 M l/7 



. 2/7 „ 

M 8 cm . 



(4) 



The typical interval between bursts is about 30 min- 
utes (Aptekar et al. 1998). We assume the bursting and 
accreting processes are in quasi-equilibrium, and the total 
energy release during a burst is 



E tot = Mx30minutcs x 30MeV 
trip 

- 2.1 x 10 46 Ms MeV 

- 3.3 x 10 40 Ms ergs. 



(5) 



The durative time during bursts ST ~ 5s (Aptekar et 
al. 1998), then the busting luminosity Lburst ~ -f^ ~ 
6.6 x 10 39 Ms ergs-s -1 ~ ^Edd- Therefore, the crust 
should be expanded very much, which could be the rea- 
son for the observed phase-lag (see section 4). The effective 
temperature T e g is 



t cS = (^r ) 1/4 



- 1.8 x 10* Bi['R e 

? l/7 p -9/28 M - 



8 R l/7 R -9/28 -1/28 ^/28 R 



(6) 



16Bi['R 6 



1/28 ■ 5/28 

' Ms keV. 
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where a = 5.67 x 10~ 5 ergs • s _1 cm _2 K -4 is the Stefan- 
Boltzmann's constant. While, in observation for GRO 
J1744-28, X-rays from 10 to 100 keV have been detected. 

Also, the column density <7 crus h of phase transition 
matter in the crush of crust should be 



^crush 



M X 30minutcs 



5.2 x 10 -B,; RI M- 



1.0 x lti*B$B?< r Mr 1/r mT 



Mq • cm 
g • cm~ 2 , 



and the extent scale h crus h of phase transition matter can 
be estimated as 

h _ £crush 

'^crusn p 

~ Q.1B{[ 7 R- 9/7 Mr 1/? tilT Pii cm«i/ crust , 

which means only a small faction of matter in the crust 
has been transformed from hadron into three-flavor quark 
matter. 

3.1. The Bursting Process is Unstable 

In equ.(3), it is said that cr max = <r m (zn = 0). In fact, as 
the temperature of the crust is not zero, ions can penetrate 
a small potential barrier with height of <~ kT, where k is 
the Boltzmann constant, T is the temperature of matter 
at the base of the crust. From equ.(A6) 



V 



3V„ 



!fV q z + 4 



(forz > 0) 



of Xu & Qiao (1998), we get the differential relation of V 
and z as 



SV 



3V„ 



(V^ z + 4 ) 2 
-4.55zn, 



'-VgSZU 



for V q = 20MeV, z u = 0. If ZSV ~ kT ~ 10~ 2 MeV 
(the charge number Z of ion is chosen to be 26), then 
Szu ~ 10~ 4 <C 1. From equ. (2), the differential relation 
of cr m and z is 

5a m ~ -8.6 x l0~ 19 5z 11 r ] p 2 1 / 1 3 A- i * /3 Z 26 MT 1 Rl 

for V q = 20MeV, z n = 0. Hence, if kT ~ 10~ 2 MeV, 
<C cr m , which means equ. (3) is a good ap- 



<5o+ 



10 



-23 



proximation for the maximum value of column density. 

When the accreting crust is heavier than that the 
Coulomb force can support, the bottom part of the crust 
should be pushed through the Coulomb barrier. The re- 
leased energy from phase transition would heat up the 
bottom crust, and according to the differential relations of 



V — z and <7 m — z, more ions should penetrate the Coulomb 
barrier. In a word, a positive feedback process happens. If 
the temperature in the feedback process (and before the 
crust have taken a global expansion) is risen to be T + 5T, 
while ST is assumed to be order of 10~ 3 MeV, the total 
mass been pushed through should be order of 10~ 24 (note: 

CTcrush ~ KT 24 ) Mq- Cm" 2 . 

4. The Phase-lag 

As discussed below, the bursting can cause matter being 
expanded large enough to drag matter to latter phase, and 
the observed phase-lag can be explained. 

The total mass of one polar cap crust M po i ar is 



M, 



polar 



4/7 



Let a be the inclination angle (Fig.l), the total rotational 
energy of the polar crust is 

Jpoiar = iM polar f7^i? 2 sin 2 a 

~ 4 - 6 v in 41 n.n 2/3 .4~ 2/3 7 M- 6/7 U 37/7 R- 4/7 

(1.2*11+4)2 X iU Wll A 56 Z 2 &M 1 R 6 ±S U 

where fl = ■jt is the angular velocity of the star. When 
bursting, the matter of crust could be expanded to a max- 
imum height r (Fig.l). Assuming this expanded crust has 
a uniform angular velocity ft, from rotational energy con- 
servation, 




Fig. 1. A demonstration of bursting process of accretion po- 
lar-crust, ft is the rotational axis, u is the magnetic axis, a is 
the inclination angle, and R is the radius of strange star. When 
bursting, the crust would expanded to a maximum height r. 



Jpoiar = g Af P oiar^ 2 (r 2 + Rr + R 2 ) sin 2 a, 



we get 



n = n 



m 2 



r 2 + Rr + R 2 ' 
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If there is a expansion in the crust, r ^ R, then Slo 7^ < 
fio- Therefore, in the corotation frame, the crust matter 
would be drag to a later phase, and the total drag kinetic 
energy is 

£k = /;i(^fdx)[(r»o-^)xsin 2 a ] 2 

= ±M polar tt 2 sin 2 a(r 2 + Rr + J R 2 )(1 - y ^jg^ ) 2 - 

While, in the corotation frame, the magnetic force for mat- 
ter to get back is F = 2 • ^ • A p since the Maxswell tense 
force is (Alfven & Falthammar 1963). So, the shifted 
distance can be estimated as 

= mh; M p°^ n ?) sin2 a ( r2 + Rr + r2 ) 

n - / „ 3K 2 \2 

~ <t££&? M sin2 a ^n 3 A^Z^M^E^B^ cm, 

(7) 

here the function /(r) is 

/(-i-p+S+fli-p-V ^Fn? )'. <«> 

The degree of phase- lag is (for zn =0), 

= 1801 x 

7T it sin a 

- 3.2° x 10 3 sina/(r)?7p 2 { 3 A 56 2/3 Z 26 Aff 1 i?|B n 2 . 

(9) 

The numerical calculation of /(r) show that could be 
consistent with observed phase-lags if r <~ 1. Oli? or r—R = 
100 meters. 

5. Conclusion and Discussion 

A model for GRO J1744-28 is suggested in this paper, with 
which many observational properties, such as the phase- 
lags, the burst energy and the characteristic photon en- 
ergy, could be explained well and simply. There might be 
two polar cap crusts of GRO J1744-28. The intermittent 
phase transition of polar crusts could result in the burst 
process when an accreting crust is heavier than that the 
Coulomb force can support. 

Bildsten & Brown (1997) argue that the bursts of GRO 
J1744-28 during the peak of the outburst are most likely 
not the thermonuclear origin. Nevertheless, if the energy 
release by phase transition is considered, the thermonu- 
clear reactions could have contribution to the bursting 
process in GRO J1744-28. Taam & Picklum (1978, 1979) 
recognized that if hydrogen burns at sufficiently high den- 
sity (~ 10 6 g cm -3 ) in the envelop of an accreting neutron 



star, helium must eventually burn within the same mass 
layer. This mixed burn of hydrogen and helium could oc- 
cur in the accreting crust of GRO J 1744-28, which might 
be triggered by the energy release of phase transition. 
Eventually, a burst appears, with a global expansion of 
^100 meters in the crust. 

Rappaport & Joss (1997) think that X-ray pulsars, 
GRO J1744-28, the rapid Burster, and type I X-ray 
bursters are in a series when the crucial parameter B 
(magnetic field) decreases. However, the magnetic field of 
the recently discovered bursting X-ray millisecond pulsars, 
SAX J1808.4-3658, is only - 10 8 Gauss (Wijnands & Klis 
1998). The spectra of GRO J1744-28 and SAX J1808.4- 
3658 are hard. The characteristic photon energy of GRO 
J 1744-28 is comparable to most other X-ray pulsars but 
higher than other X-ray bursters. On the other hand, SAX 
J1808. 4-3658 has a remarkably stable Crab-like power law 
spectrum (Gilfanov et al. 1998). Where and why is SAX 
J1808.4-3658 in the scries of Rappaport & Joss (1997)? 
We have to deal with the radiative mechanism of X-ray 
emission in GRO J1744-28 and SAX J1808.4-3658. The 
general standard pulsar high-energy emission mechanism 
could not work well. Nevertheless, a revised version of the 
standard mechanism considering the electromagnetic ef- 
fect of accretion disk could be a good candidate-model for 
the radiative process of bursting pulsars. Detailed research 
in this scenario should be necessary. 
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